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ABSTRACT 

We examine the significance of a planar arrangement in the spatial distribu- 
tion of the Milky Way's globular clusters (GCs). We find that, when separated on 
the basis of horizontal branch morphology and metallicity, the outer-most canon- 
ical young halo GC sample (at galactocentric radii in excess of 10 kpc) exhibit 
an anisotropic distribution that may be equated to a plane (24 ± 4) kpc thick 
(rms) and inclined at 8° ± 5° to the polar axis of the Milky Way disk. To quantify 
the significance of this plane we determine the fraction of times that an isotropic 
distribution replicates the observed distribution in Montc-Carlo trials. The plane 
is found to remain significant at the > 95% level outside a galactocentric radius 
of 10 kpc, inside this radius the spatial distribution is apparently isotropic. In 
contrast, the spatial distribution of the old halo sample outside 10 kpc is well 
matched by an isotropic distribution. The plane described by the outer young 
halo globular clusters is indistinguishable in orientation from that presented by 
the satellite galaxies of the Milky Way. Simulations have shown that the planar 
arrangement of satellites can arise as filaments of the surrounding large scale 
structure feed into the Milky Way's potential. We therefore propose that our re- 
sults are direct observational evidence for the accreted origin of the outer young 
halo globular cluster population. This conclusion confirms numerous lines of evi- 
dence that have similarly indicated an accreted origin for this set of clusters from 
the inferred cluster properties. 

Subject headings: globular clusters: general — Galaxy: formation 
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Introduction 



It has long been noted that the most luminous satellite galaxies of the Milky Way 



(MW) exhibit an asymmetric spatial distribution as seen on the sky (Lynden-Bell 1976 



Kunkel & Demers 1976 Lynden-Bell 1982). Kroupa et al. (2005) demonstrated that the 



spatial distribution of the 11 most luminous MW satellites is inconsistent with an isotropic 
distribution in 99.5% of random realisations. Rather, these satellite galaxies describe a 
planar structure that extends 254 kpc from the MW and possesses a scale height of ~ 20 



kpc; the MW plane of satellites (PoS) (Metz et al.|2007 2008]^ Metz et al. (2009a) extends 



this work to include the newly discovered ultra-faint dwarf galaxies with minimal change to 
the determined spatial parameters of this plane. 

The PoS so described is seen to be within ~ 8° of polar to the rotation axis of the 
MW. The Andromeda galaxy (M31) is the only other extensive satellite system for which 
we can describe the 3-d structure with meaningful accuracy. A similar degree of planarity 



is evident around M31 (Koch & Grebel 2006). Metz et al. (2009a) updated this analysis 



to include increased numbers of M31 satellites and find that the plane of satellites around 
M31 is inclined at 59° with respect to the disk (it should be noted that the number of M31 
satellites has increased by ~ 50% since the work of Metz et al. (2009b) however; see Martin 



et al. (2009), Richardson et al. (2011)). 



The configuration of satellites about their host galaxy rose to prominence with the study 



of Holmberg (1969) who proposed that satellites are preferentially found near the poles of 



the host galaxy disk (the 'Holmberg effect'). Subsequent studies based on extensive redshift 
surveys have not confirmed this effect (Sales & Lambas 2009 Bailin et al. 2008). When 



averaged over all host-satellite pairs, the satellites of relatively isolated galaxies (akin to the 
MW and M31) have a tendency to be found near the major axis of their host ( Brainerd||2005 



Agustsson & Brainerd 2010). The findings that the satellite planes of the MW and M31 are 



highly inclined are therefore unexpected. 



Recently, Kroupa et al. (2010) have argued that the existence of the MW PoS is one of a 



number of critical challenges for Lambda CDM (ACDM) cosmology. Kroupa et al. propose 



the MW satellite system is a system of tidal dwarf galaxies (TDGs) that resulted from a 
major 'wet' merger that gave rise to the galactic Bulge. However, Libeskind et al. (2010) 



and Lovell et al. (2010) contend that the observed spatial distribution and possible rotating 



configuration of the MW PoS is not an insurmountable problem for ACDM and can arise 



We do not use the term 'disk of satellites' of these authors as the term 'disk' implies a circular or 
rotationally-supported configuration in three dimensions that is not constrained in our discussion to follow. 
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naturally from hierarchical structure growth. In their simulations the MW satellite system 
properties can be frequently matched when accretion occurs preferentially along filaments 
of cosmic dark matter (DM) structure. Issues concerning cosmogony aside, the MW PoS is 
accepted to comprise a system that has resulted from accretion to the MW halo. 



It was the seminal work of Searle & Zinn (1978) that utilized the horizontal branch 



(HB) morphology of the outer halo globular clusters (GCs), together with the absence of 
an abundance gradient within the outer halo GCs, to propose that accretion onto the early 
MW was substantial, and in the case of the outer halo, the dominant formation mechanism. 



Subsequently, Zinn (1985, 1993) and more recently, Mackey & Gilmore (2004) present a 



reanalysis of the Searle & Zinn result. In these previous studies the MW GC population is 
split on the basis of HB morphology and metallicity into three sub-groups: an old halo GC 
population which have blue HB morphology at given [Fe/H] and ages ~ 13 Gyr; disk/bulge 
GCs of relatively high metallicity and exclusively red HB morphology; and a 'young halo' 
(YH) population of GCs with comparatively red HB morphology at given [Fe/H]. 



The YH GCs, as defined by Mackey & Gilmore span a distance range far in excess 
of the other populations: from 2 - 120 kpc, as distinct from, for instance, the old halo 
population where the majority reside within 30 kpc. In terms of kinematics, the YH GC 
system shows no evidence for net rotation, while the old GC system shows prograde rotation. 
Mackey fc Gilmore demonstrate the young system resembles the GC systems of the Large and 
Small Magellanic Clouds (LMC and SMC), Fornax and Sagittarius dwarf galaxies in terms 
of metallicity range, age and HB morphology. Furthermore, the YH GC population possess 
generally more extended structures than their old halo counterparts. The distribution of core 
radii for the YH system is seen to match that of the above mentioned four external dwarf 
galaxies, and distinct from the core and half-light radius distribution of the old halo GC 



subsystem to high significance (Mackey & Gilmore 2004). Based on these above similarities 



between the YH GC subsystem and the external GC systems, Mackey & Gilmore (2004) 



propose that the YH GCs are of extra-galactic origin (in accordance with the result of 



Searle & Zinn 1978), whereas the majority (but not all) of the old halo and bulge/disk GCs 



are proposed to have formed in-situ under some form of dissipative collapse as originally 



envisaged by Eggen et al. (1962). 



The study of Yoon & Lee ( |2002 ) presented evidence for a planar arrangement of a set 
of seven Galactic GCs (GGCs). These seven clusters were selected based on the RR Lyrae 
Oosterhoff dichotomy (as Oosterhoff Type-2 clusters) and metallicity (the lowest metallic- 
ity clusters in the Type-2 group). Yoon & Lee (2002) propose these clusters have been 



accreted from an external galaxy. However as pointed out by Catelan (2009), the present- 



day Milky Way satellites possess GC systems that do not match the Oosterhoff dichotomy 
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and the distribution of GGCs in the HB morphology-metaUicity plane. This discrepancy 



with present-day satelhte galaxies aside, the studies of Marin- Franch et al. (2009), Dotter 



et al. (2010), and Dotter et al. (2011) show that the age-metallicity relation of the GGCs 
possesses two branches - one with near constant old age of ~ 13 Gyrs and another branch 
to significantly younger ages with increased [Fe/H]. The younger branch of GGCs is found 
to be dominated by GCs associated with the Sagittarius and Canis Major dwarf galaxies. 



This leads Forbes & Bridges (2010) to conclude that approximately one quarter of the GGC 



system has resulted from accretion. 

In this study, we show that the YH GCs of the MW (as defined by a combination of 
age estimators) are confined to a plane that is indistinguishable from the MW PoS. The YH 
GCs act as tracers of the dwarf galaxies that have been disrupted in the formation of the 
outer halo. The consistency in spatial alignment of early accretion events to the alignment 
of current satellites suggests a preferred orientation for accretion from the local large scale 
structure that has remained consistent for a large fraction of a Hubble time. 



2. Partition of the globular cluster population by relative cluster age 



When excellent self-consistent photometry exists, such as in the study of Marm-Franch 



et al. (2009), it is possible to derive accurate relative ages even amongst clusters of advanced 



age. These authors compare the luminosity of the main-sequence turn-off to that of the 
lower main-sequence and hence express the relative age for a sample of 64 nearby GCs. 



Marin-Franch et al. finds evidence for two groups of GCs: an old GC population with an 
age dispersion of ~ 5% and no age-metallicity relation, and a young GC population with an 
age-metallicity relation similar to that defined by the GCs associated with the Sagittarius 
dwarf galaxy. 



The age determination method of Zinn (1993) relies on a relation between age, metal- 



licity, and HB morphology to resolve relative cluster age differences. |Zinn| defined a fiducial 
sequence for inner halo clusters {Rgc < 6 kpc) in a HB-typ^ against [Fe/H] diagram and 
then measured the difference between the HB-type of a given cluster and that of the fidu- 
cial sequence. Those with a difference in HB-type > —0.4 were classified by Zinn (1993) as 



young clusters, the remainder old clusters. Mackey & Gilmore (2004) has revisited the issue 



of relative age partition. Rather than using spatial information in the derivation of their 
fiducial these authors utilise the model evolutionary isochrones of Rey et al. (2001). The 



2HBR 



(B-R) 



where B is the number of stars bluer than the instabihty strip, V is the number of 



_{B+V+R)_ 

stars within the instabihty strip, and R is the number of stars to the red of the instabihty strip. 



isochrone corresponding to the canonical ancient GC population traces the those clusters 
with the bluest HB-type at a given [Fe/H] (see Figure [T|. In the study of Mackey & Gilmore 



(2004) the young halo GG population is defined as those clusters with a difference in HB- 



type in excess of —0.3 from that of the fiducial at a corresponding [Fe/H]. This corresponds 
to a minimum age difference of around -0.6 Gyr at HB-type = 0. It should be noted that 
there is an important limitation the application of this technique for estimating the relative 
ages between clusters. At the extremes of HB-type the isochrones become degenerate within 
observational uncertainties. We note that parameters other than age can also be responsible 
for the location of a cluster in these groupings. For example, in the case of NGG 2808 the 
extreme blue horizontal branch morphology is likely driven by He abundance variations. 



In the study of Yoon & Lee (2002) relative age differences of around 1 Gyr are inferred 



from an examination of the mean period of type-ab RR Lyraes {{Pab)) with cluster metallicity. 
Oosterhoff Group I clusters possess {Pab) ~ 0.55 days and are more metal-rich than those 
of Group II with {Pab) ~ 0.65 days. Yoon & Lee ( 2002[ ) demonstrate that Group II clusters 
may be further split into 'old' and 'young' clusters (Groups Il-a and Il-b, respectively). 
The benefit of this method is that we can continue our classification to the blue extreme of 
HB-type. A limitation on the method is that the {Pab) of many clusters is either limited 



by intrinsically low RR-ab numbers or the absence of studies in the literature (see Catelan 



2009). 



Our partitioning of the GGG population will rely on a combination of the latter two 
methods. In the regime where both methods are viable (that is HB-type < 0.8, and {Pab) > 
0.60 d) there are six objects in common: both methods classify five as 'young', the exception is 
the GC Rup 106 for which the {{Pab), [Fe/H]) method would suggest 'old' and the (HB-type, 
[Fe/H]) method would suggest 'young'. Examination of the MSTO in the GMD demonstrates 



that it is indeed younger than most GCs (Marm-Franch et al. 2009). For the analysis that 
follows we will use the following criteria: 



if ((HBR - EBRfiduciai < -0.3) and HBR< 0.8) 
or ([Fe/H]< -1.9 and {Pab) > 0.6 days) 
then classification = young halo cluster 

else if ([Fe/H] > -0.8) 

then classification = disk/bulge cluster 

else classification = old halo cluster 



Graphically, our partition of the GG population can be seen in Figure [1} We take the 
[Fe/H], and in the analysis to follow, the distances and positions of GGs from the compilation 
of Harris ( 1996 ) (2010 edition). The HBR is taken from the 2003 edition of the Harris catalog. 
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Koposov 1 & 2 are not included here as they lack sufficient dat to constrain their HB-type. 
The pulsation properties of the clusters are taken from the compilation of Catelan (2009). 
The three GC subgroups young halo (YH), old halo (OH), and disk/bulge (DB), as defined 
above contain 30, 56, and 28 clusters respectively. Compared to the work of MG04, our 
selection criteria sees the clusters NGC 5024, 6341, and 7099 included in our YH grouping 
on the basis of the pulsation properties of their RR Lyrae populations, whereas these clusters 
were previously assigned to the OH grouping. 



3. The spatial distribution of globular clusters 



We now investigate how isotropic the spatial distributions of the three GC populations 



are. We utilise the method described in detail in Metz et al. (2007). In brief, we utilise a 



least squares technique to define the plane that best matches the spatial distribution and 
determine the normal vector to this plane. Since the least squares method is unweighted we 
account for distance uncertainties by repeat realisations of the system with objects randomly 
shifted along their lines of sight according to a normal distribution function in which the 
variance matches distance uncertainty to each object. This provides an estimation of the 
robustness of a disk-like population. If the GC sample is not distributed in a well-defined 
plane then repeat realisations of the sample to account for distance uncertainties will show 
a large scatter in the orientation of the normal vector. A large number of realisations (10^) 
of the system are made. 



This results in a set of m unit normal vectors n for which 
components. We then define the matrix M as 



^^11 Uii ^i, 



are the cartesian 



M 



/ xi yi zi \ 

X2 m Z2 
\ i-m Vm / 



and perform an eigenvalue analysis of the matrix T=M^M. The eigenvector that cor- 
responds to the largest eigenvalue of T is the principal axis of the set of input vectors. This 
principal axis defines the orientation of spatial anisotropy in the input GC sample. There 
are three eigenvalues of T, ti < T2 < T3, where Ts is the eigenvalue of the principal axis. In 
the case of a uniform spherical distribution ti = T2 = r-^ = 1/3. Following the formalism of 



Fisher et al. (1987, their section 3.4), the following parameters derived from the eigenvalues 
provide us with a measure of how isotropic the distribution is. We define a shape parameter 
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_ In(r3/r2) 
ln(r2/ri) 



and a strength parameter C,: 



C = ln(r3/n) 



The parameters 7 and C, enable the classification of the spherical distribution of normal 
vectors. For example, a high value of 7 (say 3 or more) indicates a clustered distribution 
of normals on the sky. When combined with a high value for C,, this indicates a strongly 
condensed distribution. A low value for 7 indicates the presence of a girdled distribution, 
i.e. a distribution that forms a band around the sky. C, then describes how strongly confined 
the distribution is to this girdle. 

We now consider the distributions of our three GC groups. We compare the spherical 
distribution of the normal vectors to the best-fitting plane with those arising from sim- 
ulations of an isotropically distributed model GC population. The isotropic GC model is 
constructed from a density power law of index —3.5, which is appropriate for the overall MW 



GC population (Parmentier & Grebel 2005). In galactic coordinates this model population 
is uniformly distributed. The number of simulated GCs is matched to the observed number 
of GCs in each age group. 

The results of this comparison are seen in Figure [2j Here the shaded regions show the 
probability distribution function of the isotropic sample (realised 10^ times). The star symbol 
in each panel shows the peak of the probability density function for the observed GC sample 
(simulated 10^ times to account for the distance uncertainty to each object). In the case 
of the OH and BD GCs, the spatial distribution of the observed sample is indistinguishable 
from that of an isotropic distribution of the same size. With no restriction on the inner 
radius of the YH GC sample, there is a suggestion from the offset of the probability density 
function peaks that the YH sample is less well described by an isotropic parent distribution 
(see Figure |2] lower right panel) . 



To quantify the hypothesis that the spatial distribution of the observed sample is derived 
from an isotropic parent distribution we determine the contour lines seen in Figure [2j These 
contour lines enclose regions of the diagram which contain a given percentage of the isotropic 
sample. The percentage of times the centroid of the observed probability density function is 
replicated by the isotropic simulated sample for each GC group is given in Table [T] 
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The validity of a spatially isotropic parent population does not change significantly for 
the OH and DB GC groups as we constrain the galactocentric radius of the sample. In 
the case of the OH GC group, 19 clusters are of Rqc > 10 kpc and 6 are Rgc > 20 kpc, 
however, these samples are consistent with an isotropic parent distribution in 54% and 56% 
of trials respectively. In the case of the YH GC group, however, the hypothesis of a spatially 
isotropic parent distribution becomes increasingly unlikely as we restrict the sample to more 
distant objects. As seen in the 7-^ plane (Figure |2]) the YH GCs are increasingly located at 
larger 7 and ( with increasing galactocentric distance. This is indicative that the spherical 
distribution of the normal vectors to the best-fitting plane for the population is becoming 
increasingly confined to a preferred direction. That is to say, the YH GCs are seen to be 
increasingly confined to a planar alignment as we progress to larger galactocentric distances. 
In Table [T] we summarize our findings regarding the orientation and rms thickness of the 
plane that best matches to the observed sample. 

We now turn our attention to the spatial distribution of the satellite galaxies of the MW. 
Figure |2] shows the corresponding 7-^ plot for the 11 most luminous satellites of the Milky 



Way (i.e. identical to Metz et al. (2007)). The significance, orientation and thickness of a 
plane to the classical dwarf galaxy sample is given in Table [TJ Our results may be compared 
to those of |Metz et al] ( |2007D who find i = 158:2, b = -li:9, and a rms thickness, A=32.6 
kpc. Figure [2] (center) shows our findings with the inclusion of 10 recently discovered low- 



luminosity MW satellites (replicating the sample of Metz et al. (2009b)). Here our results 
may be compared to those of Metz et al. (2009b) who find £ = 159:7±2:3, b = — 6:8±2:3, and 



A=24.9±l.l kpc. It may be argued that on the basis of the results of Besla et al. (2010) that 



the Magellanic Clouds (MCs) be considered separately since their proper motions suggest 
they are on their first passage by the MW. On the other hand, since the MCs lie on a plane 
common to the remaining 9 'classical' satellites we propose that if the MCs are on their first 
infall that they are further evidence in favor of preferred accretion along a common plane. 
In both of the above samples we recover the significance, orientation and thickness of the 
planar distributions discussed in the literature. 

It is apparent from Table [TJthat the orientation and thickness of the plane that describes 
the spatial distribution of the YH GCs at Rqc > 10 kpc is indistinguishable from that of 
the classical MW dwarf galaxies and the recently discovered 'ultra-faint' dwarf galaxies. For 
this reason we combine these three samples to define a common plane of satellites (PoS). 

Figure [3] shows the spatial distribution on the sky of the classes of object we have 
considered. The PoS, as defined above, is shown as the solid line. The grayscale probability 
density function seen at £ ~ 156° shows the positions of the normal to this plane in 10^ 
random realisations of the input catalogue distance uncertainties. 
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Table 1: The results of our determination of best fitting plane to the distribution of objects 
as described in ^ 



GC 


Rgc 


Number 


£ 


b 


rms thickness*^ 


Significance'^ 


Group 


(kpc) 


clusters 


n 


n 


(kpc) 


(%) 


OH 


All 


56 








54 


OH 


> 20 


6 








56 


BD 


All 


28 








53 


YH 


All 


30 


144 ± 6 


-13 ±5 


28 ±5 


64 


YH 


> 10 


23 


156 ± 6 


-8 ±5 


24 ±4 


95 


YH 


> 20 


10 


169 ± 5 


-9 ±6 


23.2 ±2.7 


95 


YH 


> 30 


7 


162 ± 5 


-8 ±7 


18.5 ±2.3 


99.2 


CS'^ 


All 


11 


156.8 ±2.7 


-11.6 ±2.9 


19.3 ±2.1 


99.1 


CS+UFD'i 


All 


21 


168 ±4 


-13 ±3 


25 ±3 


98 


YL02^ 


All 


7 


183 ± 15 


-3 ±9 


19 ±4 


96 



"Calculated from repeat random simulations of the distance uncertainties for each object. See text for details. 

''Percentage of times the observed result is not replicated drawing from an isotropic population. 

'^The 11 most luminous ('classical') satellite galaxies of the Milky Way. 

''As above with the inclusion of 10 low luminosity dwarf galaxies. See text for details. 



"The cluster sample of Yoon & Lee (2002). 
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The top left panel of Figure |4] shows a section through the vicinity of the MW as seen at 
a viewing angle edge on to our derived PoS. We note that the OH GC NGC 2419 (the open 
diamond in the figure) resides well away from this plane at a distance of 83 kpc, 3.4 times 
the thickness of the PoS. There is increasing evidence that NGC 2419 is an accreted object. 



Specifically, Cohen et al. (2010), argue that due to its internal [Fe/H] spread this cluster is 
the former nucleus (or nuclear star cluster) of a disrupted dwarf. Potentially therefore NGC 
2419 had been accreted away from the plane apparently preferred by the vast majority of 
the MW's accreted objects. 



Our findings reinforce those of Yoon & Lee (2002) from a sample of GCs restricted to 



Rgc < 20 kpc. In that study, these authors isolated a set of young halo clusters based on 
metallicity and RR Lyrae pulsation period (as discussed in Their set of seven young 
GCs (their 'group Il-b' ) describe a plane orientated towards i = 183° ± 15°, b = —3° ± 9° 
(at 96% significance compared to an isotropic distribution), which is consistent with the PoS 
derived here. The PoS is therefore, persistent in the spatial distribution of the MW's YH 
GCs to the inner galaxy. 

Although we here use YH GCs as a tracer of accreted systems, it is expected that 
accreted systems will deposit both YH and OH clusters. As pointed out by |Da Costa &: 



Armandroff (1995) the Sgr dwarf is contributing to both groups. Furthermore, Mackey & 



Gilmore (2004) propose that 15-17% of the OH GCs were contributed via accretions. The 



fact that our analysis had not shown any radial anisotropy in the distribution of OH GCs is 
perhaps a consequence of the smaller fraction of accreted systems in this grouping compared 
to the YH GCs. 

We have shown that the YH GC population at galactocentric distances exceeding 10 kpc 
exhibits the same orientation (within uncertainties) as that defined by the MW satellites. 
This argues for a strong causal link between the origins of the outer halo GCs and the MW 
satellites. In the next section we discuss the origins of the planar arrangements of dwarf 
galaxies around the MW and M31, and by direct association, the origin of the outer halo 
globular clusters. 



4. The satellites of the Milky Way and M31: cosmological context 

Three scenarios are presented in the literature regarding the development of the spatial 
distribution of satellites: that the satellites are in-situ tracers of a prolate dark matter 
halo; that they are remnants from the break up of a large progenitor at early times (i.e. 
the satellites are tidally disrupted galaxies - TDGs); or they are tracers of collapse along 
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filaments of large scale structure (LSS). Cosmological simulations of dark matter halos shows 
that under ACDM cosmology halos are predominantly prolate ( Libeskind et al.|2005 Zentner 



et al. 2005), however, the degree to which the dark matter halo would have to be prolate to 



accommodate the observed anisotropy around the MW and M31 would be extremely atypical 



(Kroupa et al. 2010). Furthermore, the debris stream of the disrupting Sagittarius dwarf 



can not be modelled under the influence of such an extremely prolate dark matter halo (Law 



et al. 


2009; 


Prior et al. 


2009 


Keller et al. 


2008, 


2009 



Li & Helmi (2008) present a scenario in which all satellites fall in to the MW halo in 



1-2 groups and demonstrate that such an accretion history could account for the current 



satellite anisotropy. D'Onghia & Lake (2008) develop a similar scenario, reminiscent of the 



groups of MW satellite galaxies of Lynden-Bell (1976, 1982); Kunkel & Demers (1976), in 



which the groups formed in LMC-like dark matter halos. However, such a scenario has a 



number of problematic features as discussed in Kroupa et al. (2010), the frequency of such 
groups in the local volume is low; the group would have been comprised predominantly of 
dwarf spheroidal galaxies and this does not match to the morphology-density relation; and 



the group identity of satellites after infall is expected to be short lived ( Klimentowski et al. 



2009). 



Metz et al. (2007) present a scenario in which the satellites arise from the tidal fragmen- 



tation of a large galaxy at an early epoch. Under this scenario, the kinematics of the TDGs is 



retained from formation as a rotationally supported disk. Metz et al. (2008, 2009b) confirm 
this: the inferred angular momentum vectors for the majority of classical MW satellites are 
directed within a 30° region on the sky. 



Libeskind et al. (2005, 2010) present cosmological simulations that show that substruc- 
tures stream through LSS filaments that feed MW-sized halos. This accretion occurs from 
preferred directions on the sky and is coherent on scales from 1 Mpc to 20 kpc. From 
simulated systems that contain at least 11 luminous satellites (i.e. comparable to the MW 



system) Libeskind et al. (2010) find 30% possess satellites with orbital angular momenta 



aligned to a similar degree as the MW system. Lovell et al. (2010) find a similar result, 
namely that quasi-planar distributions of coherently rotating satellites arise naturally in the 
Aquarius simulations as subhalos fall in along the central spines of LSS filament^ 

In Figures [5] and |6] we show the orientation of the MW and M31 PoS in supergalactic 



^However we note that 



Kroupa et al. 



(2010) argues this is a selection effect citing that the MW-like 



galaxies with at least 11 luminous satellites account for 1.4% of simulated galaxies in Libeskind et al. (2010) 
and hence that only 0.4% (i.e. 30% of 1.4%) of all existing MW-like ACDM halos would host a MW-type 
galaxy with an appropriate satellite spatial distribution. 
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coordinates (see de Vaucouleurs et al. 1991 for definition). The supergalactic plane is defined 
by the distribution of galaxies ~ few Mpc about the MW. Within this volume galaxies are 
confined to a ~ 2 Mpc thick plane (seen side on in Figure [6]). Seen from above (Fig. [s]) 
the MW and M31 PoS are both highly inclined with respect to the supergalactic plane. As 
discussed in Metz et al. (2007) there is no apparent spatial connection between the MW and 
M31 PoS as the M31 PoS is inclined at ~ 55° to the MW PoS. It is, however, notable that the 
PoS of both galaxies is well aligned with an axis between the two major mass distributions 
within 20 Mpc: the Virgo and Fornax clusters ( Karachentsev et al.|2003 ). This was remarked 
upon by Navarro et al. (2004) in their study of the galaxies of the local supergalactic plan^ 
Here they found that not only are the satellites of nearby early-type galaxies preferentially 
confined to polar planes relative to their disks but that the galaxy disk is orientated along 
the major axis of surrounding LSS material. This is also the case for the MW and M31, 
where as seen in Fig. [5] and |6| the PoS are highly inclined to the stellar disks of their hosts. 



5. Discussion: the origin of outer halo globular clusters 



It is our thesis that the observed affinity between the planar distribution of outer halo 
GCs and that of the MW satellites implies a common origin. In this scenario, the outer halo 
GCs are accreted to the MW as part of galaxies from the surrounding LSS. They enter as 
entities residing in dark matter subhalos in LSS filaments that stream into the MW dark 
matter halo. Upon entry, potential subhalos may deliver the GCs from their natal dark 
matter subhalos via tidal disruption and dispersal of dwarf galaxies (akin to the simulations 



of Libeskind et al. 2010) or from the early tidal disruption of a moderate mass galaxy (like 



that envisaged by Kroupa et al. 2010). We discount the possibility of "free-fioating" GCs 



as this would imply dark matter dominance and in the case of at least one of the outer halo 



GCs (Pal 14; Jordi et al. 2009) there is no evidence for dark matter. 



We see in the spatial distribution of the outer halo GCs evidence of their accretion 
origin. This is a direct observation of the importance of accretion to galaxy formation as 
inferred from the properties of the MW GC system in the seminal work of Searle Zinn 



(1978). The only direct observations of this process of globular cluster accretion are in the 



disrupting Sagittarius dwarf (delivering at least five GCs into the MW halo: Da Costa & 
Armandroff||1995 Martinez-Delgado et al. 2002 Law fc Majewski|2010 ), perhaps in the case 



of the putative Canis Major dwarf (Martin et al. 2004) and in the globular clusters of M31 



However, we note that 



Yang et al. 



(2006); 



Agustsson & Brainerd 



(2006) find that this ahgnment does 



not appear to extend to scales less than the typical virial radius. 
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associated with tidal stellar streams in the vicinity of M31 (Mackey et al. 2010) 



Much corroborating evidence is presented for this conclusion is presented in the liter- 
ature. A particularly distinctive feature of the YH GCs is their distribution of core radii 
( Mackey fc Gilmore||2004 see also Mackey fc van den Bergh|2005 ) . The core radii (rc) of the 
YH GCs shows a very long tail to very large radii (all but one GC with > 9 pc resides in 



the Mackey & Gilmore grouping of GCs associated with the outer halo GCs). Furthermore, 
the Tc distribution of these GCs shows no statistical difference to that observed from a com- 



pilation of GCs from the LMC, Fornax and Sagittarius dwarf galaxies (Mackey & Gilmore 



2003). Together with the similarity in morphology of the horizontal branch between the 



external GCs and the YH GC grouping, this evidence leads Mackey & Gilmore (2004) to 
propose that all the YH GCs are accreted. 

The implications of the size distribution for an accreted origin of the outer halo GCs 



are elucidated in the study of Hurley & Mackey (2010). In this study, the authors use N- 



body simulations of star clusters in a tidal field to investigate the conditions required to 
produce and sustain a GC of large core radius. [Hurley Mackey] show that clusters may 
be born with a range of Tc governed by how much the cluster fills its initial tidal radius. 
When a cluster completely fills its natal tidal radius an extended GC can result (rc > 10 
pc). Conditions for a cluster to completely fill its tidal radius at birth are optimal in regions 
where background tidal forces are low: such conditions are best satisfied at large distances 
from MW-like galaxies and in dwarf systems (see discussion in Elmegreen|2008 and Da Costa 



et al. 2009). This again argues for the accretion of the outer halo GCs to the MW. 



The timing of the delivery of the outer Young Halo GCs through accretion is not con- 



strained by the physical parameters of the GCs. Gnedin (1997) describe the stability of the 



MW's GCs against two-body relaxation, tidal truncation, and tidal shocks due to passage 
through the disk and due to close proximity to the bulge. With the exception of two clusters 
(Pal 1 and Pal 13), GCs of R > 10 kpc are expected to be long-lived, with lifetimes of 
between 5 and 100 Hubble times. Hence, by this criterion, these systems could have been 
accreted into the MW at any stage over the last Hubble time. We can ask the question: 
how many satellites of a given mass are required to contribute the observed (conservative) 
22 Young Halo GCs at R > 10 kpc. Given the specific frequency of GCs (S'Ar)as a function of 
host galaxy luminosity (first introduced by Harris fc van den Bergh|[T981 , and e.g. discussed 



recently by Georgiev et al. 2010), we find this would require for example, approximately 2 
Magellanic-like {My ~ —18 with GC Sjq ~ 1) systems or 22 systems with Sculptor-like lu- 
minosities (My ~ — 11 with GC Sn ~ 70). If we then consider that the number of accreted 



clusters is likely to be supplemented by 10-12 OH clusters (Mackey & Gilmore 2004) we 



estimate that the MW may have experienced mergers with 3 Magellanic-like to 30 Sculptor- 
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like luminosity systems over the MW's lifetim^ Such deposition would contribute around 
half of the mass of the stellar halo. In order to contribute stars that are of sufficiently low 
metallicity (and sufficiently high [a/Fe]) to match those found in the halo, such accretion 
must have taken place early on in the chemical evolution of the accreted galaxies. 

6. Conclusions 

In this study we have partitioned the Milky Way's globular cluster system into groups 
according to metallicity, horizontal branch morphology, and the stellar pulsation properties of 
their RR Lyrae population. The grouping broadly align with relative ages of the constituent 
clusters. When partitioned in this way the Milky Way's young halo globular clusters (of 
galactocentric radius greater than 10 kpc) are anisotropically spatially distributed and the 
distribution can be described by a plane 24 ± 4 kpc in thickness (rms). The normal to this 
plane is directed towards {£, h) = (156° ± 6°, —8° ± 5°). In contrast the old halo and bulge- 
disk groups are consistent with an isotropic spatial distribution. The plane defined by the 
young halo population is statistically indistinguishable from that described by the satellite 
galaxies of the Milky Way {i, b = 166° ± 4°, —13° ± 3°). We argue therefore that the young 
halo globular clusters share the same origin as that of the system of satellite galaxies about 
the Milky Way. 

Simulations of the satellite galaxies of Milky Way-like systems have shown that the 
planar arrangement of satellites arises naturally from accretion of subhalos along filaments 
of the surrounding large scale structure which are directed towards the host halo. Once 
accreted, the globular clusters of the subhalos are delivered to the Milky Way halo by tidal 
disruption. 

The finding of a common spatial anisotropy amongst the Milky Way satellites and the 
young halo globular clusters is an observation of the accretion origin of the outer globular 
clusters. It is a direct confirmation of numerous lines of evidence that have shown the outer 
halo globular clusters to possess individual properties that are shared most closely with those 
found in the globular cluster systems of low mass galaxies. 

This paper arose from discussions at the Galactic Studies with the LAMOST Surveys 
workshop held at the Kavli Institute for Astronomy at Peking University, July 2010. SK and 



^When considering the number of mergers it should be noted that there is almost a factor of ten variation 
exhibited in the ^jv at given galaxy luminosity 
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Fig. 1. — HB-type versus metallicity diagram for the 114 Galactic globular clusters with 
suitable measurements in the 2003 & 2010 versions of the Harris (1996) catalogue. The 
clusters have been divided into three subsystems, as labelled, according to the criteria set 



out in the text. The overplotted isochrones are from Rey et al. (2001). The two lower 



isochrones are, respectively, 1.1 Gyr and 2.2 Gyr younger than the top isochrone. 
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Young GCs R>30kpc Young GCs R>20kpc Young GCs R>10kpc Young GCs (All) 

Fig. 2. — A (7, plot (see text for details) describing the spherical distribution of real- 
isations of the observed GC populations as indicated (star symbols, a random sample of 
100 realizations of the 10000 computed). The underlying filled color contours show the 
probability density function arising from repeat realisations of an isotropically distributed 
population. The solid lines show the locus over which (75% (red), 50%(yellow), 20%(light 
blue), 5%(dark blue), 2% and 1%) of the sample is of greater 7 and (. 
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Fig. 3. — The plane of satellites (solid line; derived from the young halo GC population 
of galactocentric distance greater than 10 kpc; see Table [l|) seen on an Aitoff projection in 
Galactic coordinates. Red, green filled circles represent the 11 'classical' luminous satellite 
galaxies of the Milky Way and the YH globular clusters [Rgc > 30 kpc) respectively. The 
open circles show the positions of outer OH clusters. The smaller blue points show the 
positions of the 'ultra-faint' dwarf satellites. Note that this figure is a two dimensional 
projection of the three dimensional distribution of satellites. For this reason the distance 
from the line (which is the plane of satellites seen at infinity) does not correspond to distance 
from the plane. The region of greyscale on the right of the figure shows the probability density 
function of the normal to this plane from repeated random realisations of the set of objects 
to account for uncertainties in their distance. The cross hair denotes the position of the 
highest likelihood position of this normal to the plane. 
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Fig. 4. — Top left: The PoS described by the young halo GC population of galactocentric 
distance greater than 10 kpc (see Table [T| as seen edge-on to the plane (solid line, surrounded 
by dashed lines showing the rms thickness, viewed from i, b = (160°, 0°)). The symbols are 
colored the same as in Figure [3] Top right: a zoom in on the inner halo, showing in green 
the YH clusters, as solid triangles the disk-bulge GCs, and as open diamonds the old halo 
GCs. Bottom left & right: the PoS seen face-on (viewed from i, b = (50°, 0°)). 
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Fig. 5. — The PoS of the Milky Way and M31 seen in the context of the supergalactic 
coordinate frame of reference. The shading of the elhpses represents orientation, darker is 
closer. In this X-Y projection the PoS of both galaxies are nearly side-on and generally 
directed along an axis joining nearby major mass concentrations in the form of the Virgo 
and Fornax galaxy clusters (15 and 12 Mpc distant from the Milky Way respectively). 
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